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We propose a method to read-out the spin-state of an electron in a quantum dot in a Voigt 
geometry magnetic field using cycling transitions induced by the AC Stark effect. We show that 
cycling transitions can be made possible by a red-detuned, circularly-polarized laser, which modifies 
the spin eigenstates and polarization selection rules via the AC Stark effect. A Floquet-Liouville 
supermatrix approach is used to calculate the time-evolution of the density matrix under the exper¬ 
imental conditions of a spin read-out operation. With an overall detection efficiency of 2.5%, the 
read-out is a single-shot measurement with a fidelity of 76.2%. 

PACS numbers: 73.21.La,42.50.-p,42.50.Hz 


I. INTRODUCTION 

Quantum information science holds great promise in 
the areas of secure communication and rapid computa¬ 
tion, but the physical components of a future quantum 
computer are still a work in progress. Any physical real¬ 
ization of a quantum bit, or qubit, requires several differ¬ 
ent single-qubit operations: initialization, manipulation, 
and read-out of its quantum stat^. There are many 
candidate systems that may act as qubits, including the 
spin degree of freedom of a single electron or hole trapped 
in a quantum dot (QD)P. Optically active transitions to 
many-body excited states allow the spin of the single par¬ 
ticle ground state to be influenced by external application 
of oscillating electric fields, such as lasers. Spin initializa¬ 
tion has been accomplished in a number of experimental 
situations involving a magnetic field in either the Fara¬ 
day configuration, where the magnetic field is aligned 
parallel to the optical axis,^^ or the Voigt configuration, 
where the magnetic field is aligned orthogonal to the op¬ 
tical axi^SHU A magnetic field in the Voigt configuration 
allows both spin initialization and coherent manipula¬ 
tion because the field modifies the polarization selection 
rules of the optical transition^2Hn]_ Statistically signifi¬ 
cant spin read-out has been achieved in both Voigt and 
Faraday configuration^^l^llSltni^ gf a cycling 

transition in the Voigt configuration makes a single-shot 
read-out of the spin-state very difficult. 

A single-shot measurement determines the state of the 
qubit faster than the back-action of the measurement dis¬ 
turbs the state. In a charged QD in a Voigt magnetic 
field, there is no optical transition that would leave the 
electron state unchanged with high fidelitji^. In contrast, 
the Faraday m agneti c field configuration results in cy¬ 
cling transitiongl^nSI^ which produce photons but leave 
the electron spin-state largely unchanged after emission. 
These cycling transitions allowed a recent demonstra¬ 
tion of single-shot spin-state read-out in the Faraday 
configuratioiP^, but the optical selection rules preclude 
arbitrary coherent spin manipulation beyond initializa¬ 
tion to an eigenstate. Therefore, in order to realize the 


three essential single-qubit operations of initialization, 
manipulation, and read-out, there is a need to combine 
the capabilities of the Voigt and Faraday configurations. 


In this Letter, we propose a scheme to read-out the 
spin-state of a single electron trapped in an optically ac¬ 
tive quantum dot using a cycling transition induced by 
the AC Stark effect of a strong optical field far detuned 
from resonance. In this scheme, a constant Voigt config¬ 
uration magnetic field allows rapid spin initialization and 
picosecond spin rotation via stimulated Raman adiabatic 
passag^^ltS^. We show below that the AC Stark effect is 
capable of modifying the allowed optical transitions re¬ 
sulting in spin-selective cycling transitions, which could 
be used for a single-shot measurement. 


Here, the AC Stark effect is induced by a circularly- 
polarized laser that is far detuned from the optical tran¬ 
sitions. In the limit of large detuning, the interaction 
between the laser and the QD does not significantly pop¬ 
ulate the excited states, b ut st ill shifts the energy levels 
coupled by that transitiorP221t221. When this AC Stark 
shift is much larger than the Zeeman splitting from the 
Voigt magnetic field, then the polarization selection rules 
are more similar to the Faraday configuration than the 
Voigt configuration. Therefore, we call the combined 
Voigt field with AC Stark shift the “pseudo-Faraday” 
configuration. Because the pseudo-Faraday configuration 
is induced by an optical field rather than a DC magnetic 
field, the system may be rapidly switched between Voigt 
and pseudo-Faraday configurations. This versatility will 
allow not only spin initialization and manipulation, but 
also single-shot spin read-out. We develop a model to de¬ 
termine the time evolution of the density matrix under 
resonant excitation in the pseudo-Faraday configuration, 
and use it to demonstrate the feasibility of a single-shot 
read-out of the electron spin state. 
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FIG. 1. (Color online) Energy level structure of a charged quantum dot with (a) no magnetic field, (b) a Faraday magnetic 
field, and (c) a Voigt magnetic field. Allowed transitions are shown with solid lines; weakly-allowed transitions are shown with 
dashed lines. 


II. DESCRIPTION OF THE SYSTEM 
A. Hamiltonian 

The ground state of the QD can be the empty state, 
where the lowest lying conduction and highest lying va¬ 
lence bands are empty of electrons and holes respectively. 
Or the ground state can contain a single charge, either 
an electron or hole. If the QD is in a diode structure, 
the QD charge can be stabilized and adjusted depend¬ 
ing on the bias voltage applied to the diode contact^. 
In the case of an n-i-Schottky diode structure, the QD 
can have a single electron trapped in the bound conduc¬ 
tion band state. It is the spin-state of this trapped elec¬ 
tron that may serve as a qubilP. The single-electron z- 
projection spin states are optically coupled to charged 
exciton (trion) spin states comprising a pair of electrons 
in the conduction band and a single heavy-hole in the va¬ 
lence bancP^. Because the electrons form a singlet state, 
the spin of the trion is determined solely by the spin of 
the hole. One of the electrons may recombine with the 
hole, emitting a photon and returning the QD to the sin¬ 
gle electron ground state. We name the relevant eigen¬ 
states of the negatively charged QD as \e,z+), |e, 2—), 
\t,z+), \t,z—) where e means the single electron state, 
t means the trion state, and is for the z-projection 
of the spin. Due to conservation of angular momentum, 
each trion state has an allowed transition only to the 
electron state of matching spin, and the transition to the 
opposite-spin electron state is only weakly allowed due 
to slight light-hole/heavy-hole mixin^^^^IMl xhe energy 
level structure of the charged QD is shown schematically 
in Fig. [ 3 a). 

The Hamiltonian of a negatively-charged QD in both 
a magnetic field and an electric held is: 

H = IiaJo(o’+o’+ + crlcr_) — fl ■ B — d ■ E (1) 

where cug is the transition frequency, (t_|_ and (t_ are the 
lowering operators for the z-projection spin-up and spin- 
down manifolds, respectively, fl and d are the magnetic 
and electric dipole operators, and B and E are the mag¬ 
netic and electric held amplitudes. The lowering opera¬ 


tors are dehned in terms of the electron and trion spin- 
states as 

a+ = \e,z+) {t,z+\ , , 

cr_ = ie,z-) {t,z-\ 

We decompose Ed into “atomic” {Ha), magnetic dipole 
{Hz), and electric dipole {Hjj) components {H = Ha + 
Hz + Hd)- With no magnetic held, the z± spin projec¬ 
tion states are degenerate. A Voigt conhguration mag¬ 
netic held is perpendicular to the propagation direction of 
the emitted light, which is typically the z-direction, nor¬ 
mal to the sample surface. In that case B = xB^, and 
the magnetic dipole, or Zeeman, Hamiltonian become^^ 

Hz = d‘BBx{ge,x{s\ + Se) + gh,x{s\ + Sh)) (3) 

where /is is the Bohr magneton, g^^x and gh,x are the 
electron and hole g-factors for a magnetic held in the x- 
direction, and and Sh are the electron and hole spin-hip 
operators, respectively. The spin-hip operators couple 
states of similar charge conhguration but opposite spin. 
In terms of the electron and trion spin-states, the spin- 
hip operators are dehned as 

Se = |e,z-) (e,z-F| , , 

Sh = \t,z-) {t,z+\ ^ 

The form of the Zeeman Hamiltonian in the Voigt con¬ 
hguration leads to eigenstates that are superpositions of 
the zero-held spin-states. 

To describe the AC Stark shift, we assume an oscilla¬ 
tory form for the electric held, as in a single-frequency 
laser beam, and perform the standard rotating wave 
approximatiorP^ to obtain a Hamiltonian for the unper¬ 
turbed QD and the electric dipole interaction: 

Ha + Hd =-{/iAi(tT^cr+ — a+a\_ + ctIct- — 

+d{Ei+a+ + -f H.c.)} ( 5 ) 

where Ai = uiq — uji is the detuning of the laser frequency 
wi from the QD resonance wq; d is the dipole moment of 
the transitions; 61+ and Ei_ are the complex amplitudes 
of the left and right circularly polarized components of 
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the electric field; and H.c. means Hermitian conjugate. 
We can further simplify this expression by introducing 
the complex Rabi frequencies associated with the two 
circularly polarized components: ^1+ = dEi+/h, and 
ill- = dEi-/h. We are primarily interested in circularly- 
polarized light because it will produce an AC Stark shift 
that reduces the effect of the magnetic coupling between 
the spin-states, as we will demonstrate below. There¬ 
fore, we henceforth assume left-circularly polarized light 
(rii- = 0), and without loss of generality we can treat 
rin- as purely real. 

Using Hji, Elz and Hjj from above, we can now write 
the Hamiltonian from Eqn. Q of the charged QD system 
in the rotating frame as 

Hq =-/iAi(tT^cr_|_ — cr+cr^ -|- cr^cr_ — (T_cr^l) 

-\-^BBx{ge,x{s\ A A 9h,x{s\^ + Sh)) 

+ a+) ( 6 ) 

In the basis of the unperturbed QD eigenstates in the 
rotating frame, the matrix representation of the Hamil¬ 
tonian is 


■-/iAi/2 nsBxge.x hV,i +/2 0 

^ gBBxge,x -hAi /2 0 0 

° SH 1+/2 0 SAi/2 -gBBxgh,x 

0 0 -EBBxgh,x hAi /2 _ 

(7) 

and the unperturbed eigenstates are represented as vec¬ 
tors 
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( 8 ) 


The interaction represented by and Ai is that of the 
far detuned laser that will cause the AC Stark effect. 


B. Zeeman Effect 

In the presence of a magnetic field in the Voigt geom¬ 
etry with Bx > 0 (and no laser field), the eigenstates 
are no longer the z± projections of the spin, but the x± 
projections, which are superpositions of the z-projection 
states. We can diagonalize the system Hamiltonian Hq 
with Hi_|_ = 0 and Ai = 0 to obtain the eigenvalues 

“ EBBxge.x 
^2 —EBBxge,x 
A3 = — EBBxgh,x 
A4 —f^BBxgh,x 


and eigenstates 



■ 1 ■ 


T 

S’a;-) = ^ 

-1 

0 

A 

II 

1 

0 


. 0 . 


0 . 



'O' 


■ 0 ■ 

- V2 

0 

1 


0 

1 


.1. 


-1. 


These representations of the eigenstates are in the z± 
basis, thus we can see that the xE projections are super¬ 
positions of the z-projection states. The left side of Fig. 
shows the Zeeman splitting of the electron and trion en¬ 
ergy levels for a Voigt geometry field. In a Voigt config¬ 
uration, transitions from either trion spin-state to either 
electron spin-state are allowed, as depicted schematically 
in Fig. [^c). There are no cycling transitions that might 
allow a single-shot fluorescence measurement of the elec¬ 
tron eigenstate. 


C. AC Stark Effect 

In the absence of a magnetic field and in the large 
detuning limit, the energy levels coupled by the elec¬ 
tric dipole interaction are modified by the AC Stark 
shiffPSHSU^ We can see this by determining the eigenvalues 
of the Hamiltonian Hq from Eqn. Q with Bx = 0 : 

Ai 

^ ^ A 

A2 —-Ai 
A3 “2^^ 

A ^ A 

A4-2A1 

where Wi = ^A\e is the generalized Rabi fre¬ 

quency for the far-detuned laser. The eigenvectors of the 
Hamiltonian, still expressed in the basis of the unper¬ 
turbed states, are 
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The eigenstates corresponding to the z± manifold are 
no longer purely electronic or trionic, but a superposi¬ 
tion of both. This occurs because they are coupled by 
the (T_|_ polarization of the laser. If we make the assump¬ 
tion that the detuning Ai is much larger than the Rabi 
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FIG. 2. (Color online) Splitting of the electron and trion spin 
states due to the Zeeman effect (left) and the AC Stark effect 
(right). The trion eigenfrequencies are plotted relative to the 
zero-field transition resonance frequency ojo/Stt (dashed), for 
a fixed laser detuning of Ai/27r = 2000 GHz. 


frequency fii-i-, then to first order in the ratio f2i+/Ai 
the generalized Rabi frequency is 

w.»a. + M± 0) 

and we can approximate the eigenvalues and eigenvectors 
as follows. 
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( 10 ) 


Note that these eigenvectors are normalized only to 
first order in Oi_|_/Ai. In the large detuning approxima¬ 
tion, one state, vi, is more electron-like and the other, ua, 
is more trion-like. Because the Hamiltonians in Eqns. ([^, 
and Q are expressed in the rotating frame in order 
to make the rotating wave approximation, the eigenvec¬ 
tors in Eqn. ( | 10 [ ) are also expressed in the rotating frame. 
Therefore, to determine the energies of the states we must 
add h{u}o -I- a;i)/2 to the eigenvalues for the trion-like 
states, Aa and A4, and add h{uJo — 0Ji)f 2 to those for the 
electron-like states, Ai and A2: 
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( 11 ) 


We can see from the above expressions that a circularly- 
polarized far-detuned laser results in a spin-selective AC 
Stark shift from the unperturbed energies. For the case 


of a a+ polarized field that we consider here, the z-\- man¬ 
ifold states are shifted {Ei & Ea), while the z— manifold 
states are not {E^ & £"4). For red-detuning (Ai > 0 ), 
the electron z+ energy shifts downward and the trion 
z+ energy shifts upward by the same amount. For blue¬ 
detuning (Ai < 0), the energy shifts would be the oppo¬ 
site. 

The purpose of the AC Stark laser is to shift the en¬ 
ergies without populating the trion-like states. A red- 
detuned laser is preferred over a blue-detuned one be¬ 
cause of the lower probability of inelastic absorption for 
a red-detuned laser. Inelastic absorption is the absorp¬ 
tion of a photon combined with either emission or ab¬ 
sorption of a phonon from the crystal lattice. Inelastic 
absorption of a red-detuned photon would require the si¬ 
multaneous absorption of a phonon, which is improbable 
due to the low temperature at which optical experiments 
on QDs are usually performed. Inelastic absorption of 
a blue-detuned photon, however, would require only the 
emission of a phonon, and that process can still occur 
even at zero temperature. Therefore, a red-detuned AC 
Stark laser causes the energy states to shift while mini¬ 
mizing the probability of exciting the QD. Subsequently 
we will assume that the AC Stark laser is red-detuned. 

The amount of state mixing in the z+ manifold is pro¬ 
portional to r2i+/Ai, and can thus be reduced arbitrar¬ 
ily by increasing the detuning. The energy shift, how¬ 
ever is proportional to ftf_f_/Ai. Thus, with enough laser 
power we can have a situation where the state mixing 
can be made negligible while maintaining a non-zero en¬ 
ergy shift. The right side of Fig. shows the shifting of 
the energy levels as a function of Rabi frequency at 
a fixed detuning of Ai/ 27 r = 2000 GHz. At fixed large 
detuning, the energy shift is quadratic in and for 

CT4- polarization it only affects the z+ manifold states. 


D. Pseudo-Faraday Configuration 

The Voigt geometry magnetic field couples the z+ spin 
states to the z— spin states, while the AC Stark effect 
of a (7+ polarized laser shifts only the z+ spin states. 
In the presence of both a magnetic field and a strong, 
far red-detuned ct+ polarized laser, the spin projection 
of the eigenstates and the polarization selection rules re¬ 
sult from a competition between the magnetic coupling 
between the z+ and z— states and the AC Stark shift of 
the z+ states. When the AC Stark shift is very large com¬ 
pared to the magnetic coupling, then the magnetic field 
has little effect on the system. We call this situation the 
pseudo-Faraday configuration because the energy struc¬ 
ture, eigenstates, and polarization selection rules are sim¬ 
ilar to the Faraday magnetic field configuration where 
B = zBz- We can calculate the energy levels from the 
eigenvalues of the Hamiltonian Hq from Eqn. (0, but 
because that representation is in the rotating frame, to 
obtain the eigenenergies we must perform the same op¬ 
eration as for Eqns. above. Because the pseudo- 
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FIG. 3. (Color online) Evolution of eigenstate frequen¬ 
cies as a function of first magnetic field, and then Rabi fre¬ 
quency with a fixed magnetic field of 0.1 T and detuning 
Ai/27r = 2000 GHz. The diagram to the right shows the 
energy level structure in the pseudo-Faraday configuration. 
Allowed transitions are shown as solid lines and labeled with 
their polarizations. Weakly-allowed transitions are shown as 
dashed lines. 


Faraday eigenvalues are roots of the 4 th-order character¬ 
istic polynomial of Hq they are very complicated, thus 
we forgo including them here. Instead, we calculate them 
numerically and discuss specific points of interest. 

Figure depicts the evolution of the QD eigenstate 
frequencies (oc energy) for finite magnetic field and large 
detuning, Ai ^ ^1+. In the first section (left of the ver¬ 
tical line) r2i+ = 0 and the magnetic field increases, caus¬ 
ing the previously degenerate electron and trion levels to 
split and the eigenstates to become the x±. spin projec¬ 
tions, which are superpositions of the zi spin states. In 
this configuration, transitions from either trion spin-state 
to either electron spin-state are allowed. In the second 
section (right of the vertical line), the magnetic field is 
fixed at O.I T and is increased. The AC Stark ef¬ 
fect shifts the z+ components of the electron and trion 
states, but not the z— components. When the AC Stark 
shift becomes much larger than the Zeeman splitting, the 
eigenstates are more like those of the Faraday configura¬ 
tion, except that the z— manifold is unperturbed. The 
final energy-level configuration is shown to the right of 
the plot; the allowed transitions are circularly-polarized 
and shown with solid lines, while weakly-allowed transi¬ 
tions are shown with dotted lines. 

The true Faraday configuration has cycling transitions 
that allow a single-shot measurement of the electron spin- 
stat^i^. As depicted in Fig. [^b), the two trion states of 
the QD each have two possible transitions, and the ra¬ 
tio of their probabilities - or emission rates - is called 
the branching ratio. In a cycling transition, the branch¬ 
ing ratio is very different from unity, meaning that the 
spin-preserving transition is far more likely than the spin- 
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FIG. 4. (Color online) Evolution of the branching ratio 
from the Voigt configuration = 0.1 T) to the pseudo- 
Faraday configuration. In the pseudo-Faraday configuration, 
the branching ratio is 0.02. 


flipping transition. The branching ratio quantifies the 
measurement back-action of the cycling transition and 
will determine the number of cycles that can be used for 
a single-shot measurement. 

To obtain the branching ratio, we take the ratio of 
the spin-flipping transition rate to the spin-preserving 
transition rate for one of the trion-like states: 


KV^l |d|V:4)|^ 

|(V'2 |d | V'4)|^ 


( 12 ) 


where |^i) are the pseudo-Faraday eigenstates that are 
the eigenvectors of the Hamiltonian Hq in Eqn. ([^, and 
d is the dipole moment operator (see the appendix for 
the derivation of Eqn. (HI)- 

In Fig. 1 ^ we plot the branching ratio as a function 
of Rabi frequency for fixed magnetic field and typical 
g-factorJ^ 5 e,a; = 0 . 47 , gh,x = 0 . 24 . As the Rabi 
frequency increases, the branching ratio reduces signifi¬ 
cantly, becoming similar to that of the Faraday configura¬ 
tion. For a magnetic field of O.I T, detuning of Ai/ 27 r = 2 
THz, and Rabi frequency of Di_|./ 27 r = 200 GHz, we cal¬ 
culate the branching ratio to be 0 . 02 . Measured values of 
branching ratio for solid-state systems range from O.OOl^^ 
to 0.0#SI. Our prediction of the branching ratio for a QD 
in the pseudo-Faraday configuration is within this range, 
suggesting that it may be possible to perform a single¬ 
shot read-out of the electron spin-state via resonantly 
excited fluorescence. 


III. SPIN READ-OUT OPERATION 
A. Overview 

We have demonstrated that applying a strong, far red- 
detuned, circularly-polarized laser to a charged QD in 
a Voigt configuration magnetic field results in a situa¬ 
tion similar to a Faraday magnetic field. We now discuss 
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how a practical spin read-out scheme would work in the 
pseudo-Faraday configuration. To perform a spin read¬ 
out operation, we need to excite the system in a spin- 
selective manner and detect the fluorescence. Detection 
of a photon would correspond to the electron being in a 
certain spin state, and the fidelity of the measurement 
depends on the branching ratio and the spin-selectivity 
of the excitation. 

In the pseudo-Faraday configuration, there are two 
spin-preserving transitions that are non-degenerate and 
two weakly-allowed spin-flipping transitions that are 
nearly degenerate (see Fig.[^. The spin-preserving tran¬ 
sitions are circularly-polarized, (t+ and (T_, and can both 
be excited by linearly-polarized light. The energy differ¬ 
ence between the two spin-preserving transitions, how¬ 
ever, makes resonant excitation of one of them a spin- 
selective excitation even when linear polarization is used. 
Thus, spin-selective excitation may be accomplished by 
a linearly-polarized laser tuned to resonance with one 
of the spin-preserving transitions. The fluorescence will 
be circularly-polarized and may, therefore, be distin¬ 
guished from the linearly-polarized laser scattering by 
cross-polarized detectiorpl. Alternatively, we can use a 
modal discrimination method, as in Refs. |36j and 
wherein the resonant laser is introduced into the waveg¬ 
uide mode of a planar microcavity that confines the laser 
scattering, while the QD fluorescence is emitted into the 
orthogonal cavity mode. 

After the electron has undergone some operations (e.g. 
initialization and/or manipulation) in the Voigt config¬ 
uration, the AC Stark laser can be applied to trans¬ 
form to the pseudo-Faraday configuration where read¬ 
out will occur. The application of the AC Stark laser 
can be rapid compared to the switching time of a mag¬ 
netic field, but if the laser field is applied too fast then 
the transition from Voigt to pseudo-Faraday will be non- 
adiabatic. Prior to the application of the AC Stark laser, 
the electron may be in any arbitrary superposition of 
the two Voigt eigenstates, which are the x± projection 
states of the spin. The transition from Voigt to pseudo- 
Faraday configurations will occur adiabatically if the AC 
Stark laser is turned on slowly relative to fi/6e, where 
6e = 2fiBge,xBx is the energy splitting of the electron 
spin states in the Voigt configuration. For a 0.1 T mag¬ 
netic field and typical electron g-factoil^ ge.s = 0.24, the 
electron spin precession period is h/6e = 120 ps. There¬ 
fore, the AC Stark laser must have a rise time of about 
1 ns or greater. If the adiabatic condition is satisfied, 
the population of the xzt state in the Voigt configuration 
transitions without change to the population of the 
state in the pseudo-Faraday configuration. This is for 
the case of a cr+ polarized AC Stark laser; for a cr_ po¬ 
larized laser, the mapping from x-basis to z-basis would 
be the opposite. Due to the adiabatic mapping of the 
Voigt x-basis to the pseudo-Faraday eigenbasis, detection 
of the spin-selectively excited fluorescence in the pseudo- 
Faraday configuration is equivalent to a projective mea¬ 
surement of the spin onto the x± states in the Voigt con¬ 


figuration. Furthermore, a measurement in any arbitrary 
spin projection basis can be performed by preceding the 
AC Stark laser with a rotation of the electron spin Bloch 
sphere via stimul ated Raman adiabatic passage by addi¬ 
tional laser pulsed^lUJ The capability of using multiple 
measurement bases for identically prepared states would 
allow full quantum state tomographySSl to be performed 
on the electron spin state. 


B. Floquet-Liouville Supermatrix Approach 

To demonstrate the feasibility of a single-shot read¬ 
out of the electron spin state we numerically calculate 
the time evolution of the QD density matrix p{t) under 
the conditions outlined in the previous section. This re¬ 
quires the addition of another electric dipole Hamiltonian 
describing a second, near-resonant laser at a frequency 
UJ 2 that can spin-selectively excite population from the 
electron states to the trion states. In the rotating frame 
with rotation frequency wi, the Hamiltonian of the sec¬ 
ond laser interaction is oscillatory: 

■ 0 0 0 0 - 

_ ft 0 0 0 0 

2 ^ 2 + 0 0 0 

.0 0 0 . 

where A 2 = ujq —102 is the detuning of the near-resonant 
laser, and D 2 + = dE 2 +/h and ^ 2 - = dE 2 -/h are the 
complex Rabi frequencies for the two circularly polarized 
components of the electric field. The total Hamiltonian of 
the charged QD system, including magnetic and electric 
field interactions, can be expressed as 

H{t) = -f ^(-l)g-i(A2-Ai)i (-^4^ 

where Ed^-^'^ is the time-independent EIq from E qn. Q 
and are the two constant matrices in Eqn. ( |l3[ ). 

The time evolution of the density matrix can be de¬ 
termined by solving the Liouville equation, which can be 
extended to include spontaneous transitions by using a 
Lindblad superoperator 

|p(t) = -^[H(t),p(t)]+£(p) (15) 

In solving the Liouville equation numerically using a 
standard differential equation solving algorithm, the os¬ 
cillatory nature of the Hamiltonian H(t) requires the 
integration time-step to be much smaller than the os¬ 
cillation period 27r/(A2 — Ai). The interesting system 
dynamics, however, occur on a time scale much longer 
than the oscillation period. Thus it is inefficient to 
use this form of the Liouville equation to numerically 
solve for the long-term dynamics of the system. Instead, 
we solve the equation using a Floquet-Liouville super¬ 
matrix approaclP^, which we describe here for the spe¬ 
cific case of interest. Similar Floquet theory approaches 
have been used to describe the spectrum of resonance 


e*(A2-Ai)‘Tij.c. (13) 
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fluorescence from QDs under bichromatic near-resonant 
excitatiorPi®!!. Here, however, we need not the emis¬ 
sion spectrum but the density matrix evolution. The 
Floquet-Liouville supermatrix approach allows analyti¬ 
cal solutions to any order of approximation, meaning 
that the density matrix at any time can be calculated 
without needing to calculate all the intervening density 
matrix values. 

The unperturbed electron and trion eigenstates of Hq 
form a complete orthonormal basis {|a)} in the Hilbert 
space of the QD states. In this basis, the density ma¬ 
trix operator can be expressed using the matrix elements 
Pafiit) = {o:\p{t)\l3): 


This expression is often considered a matrix in the state 
basis {|a)}, but it can also be considered a supervector 
p{t) in the operator basis {|a) (/3|} with elements papit). 
This interpretation unfolds the 4-by-4 density matrix into 
a 16-dimensional density supervector. A similar trans¬ 
formation can be performed on the Liouville equation, 
which is a 4-by-4 matrix equation in the Hilbert space of 
the QD states. Utilizing the completeness of the {|Qf)} 
basis, 

l = ^|a)(a| (17) 


P{t) = '^Pcpit) |a) (/3| 

a/3 


(16) we can express Eqn. (15) as 


'y PaP leu) (/3| y ) I , y ) {.HakPkp PakHkp) + (■^(p))ct/3 j 1*^) ' 


at ^ \ h 

0.0 o0 \ k 


(18) 


For clarity we have dropped the explicit time depen¬ 
dence of H{t) and p{t). While the Hamiltonian can be 
expressed as a 4-by-4 matrix operating on the density 
matrix, as in the above equation, the Lindblad super¬ 
operator cannot. C{p) represents population relaxation 
processes (e.g., spontaneous emission) and decoherence 
processes (e.g., pure dephasing and all population relax¬ 
ation). Population relaxation, or Ti processes, are the 
a = P elements of {£{p))ap' 


(U(p))Q,a — y ) ( £aqPaot T £qaPqq') (^^) 

9 


where F^^ is the spontaneous transition rate from state 
|a) to state \P). The first term in Eqn. (19) represents 
transitions from state |a) to all the other states; the sec¬ 
ond term is transitions to state ja) from all the other 
states. The population relaxation rates form a 4-by-4 
matrix with elements T^p that is not generally symmet¬ 
ric: transitions from trion states to electron states occur 
spontaneously, but not the reverse. Examining Eqn. (19), 
the effects of the diagonal elements Foo, cancel out, and 
therefore without loss of generality we can set them all 
to zero. Decoherence, or T 2 processes, are the a P 
elements of {C{p))ap'- 



E 


(flaij + r^sq-) 



PaP 


( 20 ) 


where ')ap is the pure dephasing rate for the coherence 
Pap between the states jo:) and \P). The terms in the 
sum are the decoherence caused by population relax¬ 
ation. The second term is pure dephasing, also called ho¬ 
mogeneous broadening. The pure dephasing rates form 
a 4-by-4 matrix with elements ^ap that is symmetric: 

lap — 1 Pa¬ 
in contrast with the Hamiltonians, the Lindblad su¬ 
peroperator cannot be expressed as a single 4-by-4 ma¬ 
trix operating on the density matrix. However, with the 
above decompositions of L{p) and identities such as 

PaP — y ) PuP^v (^^) 

flU 

where Sj-j is the Kronecker delta, we can rearrange 
Eqn. ( |18[ ) into a supermatrix equation that allows the 
both the commutator with H (t) and the Lindblad super¬ 
operator to be expressed as a single 16-by-16 supermatrix 
operating on p(t): 


lp(t) = -^-Lit)pit) ( 22 ) 

Or, expressed using the supervector and supermatrix el¬ 
ements: 


m 


Papip) 


i 

h 


y ^ Lap-,fivit)Pfiu{t) 


( 23 ) 


{■^{p))aP 


The elements of the Liouville supermatrix L(t) are: 







( 24 ) 


(oi — /3) -^iya(0^a/i) “t" ^^ao^aa^aL' 

9 ; 


(a ^ /3) Lcp-^^{t) = {Haf,[t)5i3v - H^i3{t)Sa)j.) +ih \ X) (^ag + ^fSq) “ 7a/3 <^o 


Although in Eqn. (22) the Liouville equation is now 


expressed as an ordinary differential equation with a 
single matrix, the matrix still has an oscillatory time- 
dependence. Similar to how the Hamiltonian is separated 


in Eqn. (14), we can separate L(t) into a constant term 
and two oscillatory terms: 

L{t) = -t L(i)edA2-Ai)t ^(-i)g-i(A2-Ai)t ('25) 


where contain the constant part of H(t) and all of the 
relaxation terms from £{p), and contains only the 

Hamiltonians In detail, the supermatrix terms in 

L{t) are: 


(a = P) 

(a 7^ P) 


r(±l) _ 

Ojjv Oafi 


rw> _ ttK'JIs: _ ttK'J)/: 


(0) 

-( 0 ) 


^(0)^ 

r(o); 




- UCt t 

r(o); 




(26) 


2 7ck/3 1 

9 / 


The Hamiltonian and thus the Liouville supermatrix 
both oscillate at a frequency v = A2 —Ai. Therefore, the 
density supervector p(t) will have oscillatory components 
at frequencies that are harmonics of u. We can use a 
Floquet expansion to express the supervector as a sum 
of slowly-varying supervector coefficients multiplied by 
oscillatory functions: 


We also define operators on the Fourier space: 


+00 

Fz \m) =m\m) = J2 1*^) (*^1 

n— — oci 
+00 

Fm \n) = \n + m) Fm = \n + m) {n\ 

n— — oo 


(30) 


+ 00 


P{t) = ^ P 




We* 


(27) 


m— — oo 


Substituting Eqns. (25) and ([2^ into Eqn. (22) results in 


an infinite series of coupled linear differential equations 
for the supervector coefficients: 


E 


dt 


F irrny/P"^\t)^ e* 




pdp) „-i{n.+p)t 


(28) 


We invoke single-mode Floquet theorj® to simplify 
this infinite series of equations. To describe the rapidly 
oscillating factors, we define a Fourier state space = 
{|oo) ,..., |1) , |0) , | —1) ,..., |— 00 )} where the state |m) 
represents oscillation at the m**' harmonic of the Hamil¬ 
tonian oscillation frequency v: 


{t\m) = e*' 


(29) 


The supermatrix L(t) is an operator on the 16- 
dimensional Hilbert space defined by B(^h) = {!«) 

We define the Floquet space as the tensor product be¬ 
tween the Hilbert space and Fourier space: 0 

B(^h)- Using Floquet space, we can express the finite¬ 
dimensional time-dependent Liouville supermatrix evo¬ 


lution equation (22) as an infinite-dimensional but time- 


independent Floquet-Liouville supermatrix equation: 


d . i / X 

= -^LFPFyt) 


where the Floquet-Liouville supermatrix is 


(31) 


+00 


Lp — 


^ (Fr,®L^^^)+hy{Fz®I(H)) (32) 


and /(jj) is the identity operator in Hilbert space. Ex¬ 
pressed as a matrix in Fourier space, Lp is 
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Lp 


0 

■. L(o) + 2vIi^H) 

0 + vI(H) 

0 0 

0 0 0 

0 0 0 

0 0 0 


Each element of the above expression is a 16-by-16 Liou- 
ville supermatrix in the operator basis of Hilbert space. 
The Floquet space density supervector is an infinite¬ 
dimensional vector: 

Mt) = \m) = \m)(^\a) (/3| (34) 

m moL^ 


Expressed as a vector in Fourier space, pF{t) is 


Ppit) 




(35) 


Because the Floquet-Liouville supermatrix in Eqns. (32) 
and (33) is time-independent, the solution to Eqn. (31) 


is simple and well-known: 


pF{t) = e 


_ -iLpt/h 


Pf{0) 


(36) 


OF{t) =[■■■ 


Finally, the Hilbert space density supervector at an ar¬ 
bitrary time with arbitrary initial conditions can be de¬ 
termined from the expression 

m = OFmF{t)m<s>m) ( 40 ) 


Though the supermatrices involved are in principle 
infinite-dimensional, in practice the Fourier space only 
needs a few dimensions, after which the space can be 
truncated without significantly altering the results of the 
computation. The solution p{t) in Eqn. (40) can be cal¬ 
culated in a computationally efficient manner because it 
requires no numerical integration, only specification of 
the initial conditions and matrix algebra. 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

L(o) 

L(i) 

0 

0 

L(-i) 

1 

o' 

L(i) 

0 

0 

L(-i) 




0 0 


where pf{0) is the initial Floquet supervector, which can 
be expressed in terms of the initial Liouville supervector, 
p(0), as 


Pf(0) = (|0)®p(0)) 


0 

m 

0 


(37) 


We can now determine the Floquet space density su¬ 
pervector at any arbitrary time simply by applying the 
evolution operator Upit) = to the initial condi¬ 

tions pf{0)- In the end, however, what is needed is the 
Hilbert space density supervector, which can be obtained 
from pF{t) by applying the operator 

= (38) 

n 

The operator Opit) traces over just the Fourier space, 
leaving the Hilbert space density supervector. It is an 
infinite row vector in Fourier space, whose elements are 
Hilbert space supermatrices: 


r, . T, ^ ^ — r, ^ 

^(H) e l(^H) e l(^H) 


(39) 


C. Simulations 


We can now numerically calculate the density matrix 
evolution under experimental conditions. A spin read¬ 
out operation will consist of resonant excitation in the 
pseudo-Faraday conhguration where the QD is in one 
of the two electron-like eigenstates, z+ or z—. Each 
electron-like eigenstate has one strong and one weakly- 
allowed optical transition. To perform a projective mea¬ 
surement on the z— eigenstate, for example, the resonant 
laser is tuned to the allowed a- polarized transition (see 
Fig.i and the fluorescence is detected. Detection of 
a photon is interpreted as confirmation that the electron 
was in the z— eigenstate. The fidelity of the spin read-out 
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depends on the specificity of the excitation and fluores¬ 
cence. If the resonant laser is tuned to the tT_ polarized 
transition, but the electron is in the z+ spin state, pho¬ 
ton emission may still be detected, giving an erroneous 
signal. 

To predict the fidelity of the spin read-out opera¬ 
tion, we simulate it using typical QD parameters and 
a linearly-polarized resonant laser tuned to the cr_ cy¬ 
cling transition. The experimentally controllable pa¬ 
rameters used are: = 0.1 T, f7i+/27r = 200 GHz, 

Ai/27r = 2 THz, Q2+/27r = = 0.5 GHz. 

Note that the AC Stark mixing parameter is kept small 
(f2i+/2Ai = 0.05), but the AC Stark shift is still large: 
r2f_,_/2Ai = 27r(10GHz). In order to have the excitation 
laser in resonance with the cr_ cycling transition, the 
detuning A2 is chosen to match the difference between 
the eigenvalues of the time-independent pseudo-Faraday 
Hamiltonian that correspond to the z— electron¬ 

like and trion-like states. The QD parameters would 
not be controllable in practice, except by the choice of 
QD, but typical values are used here. The electron and 
hole g-factor^ are gf.^^ = 0.24, = 0.47. The popu¬ 

lation relaxation rates used in the simulation are 
all zero except for: the spontaneous emission rate^^ 
Fai = r42 = 1.54 GHz; the weakly allowed emission 
rateJi^r4i = r32 = 3.42 MHz; and the electron spin de¬ 
cay rate^ r2i = ri2 = 50 Hz. The dephasing rates 
used are: trion dephasin^^yai = 713 = 742 = 724 = 1.72 
GHz; electron spin dephasin^^ 712 = 721 = 12.6 MHz. 

We numerically calculate the density matrix evolution 
for two initial conditions: the z— and z-\- electron-like 
states. The photon emission rate is proportional to the 
population in the trion states and is given by: 

R{t) = r3iP33(t) -I- r42P44(t) (41) 

The average number of detected photons is the overall 
detection efficiency e multiplied by the time integral of 
R{t) over the duration T of the detection window: 

D{T) = e [ R{t)dt (42) 

Jo 

Figure [^a) shows the photon emission rate R{t) as a 
function of time after resonant excitation begins. When 
the initial state is the z— electron state (solid red curve) 
the emission starts strong and decays as the population is 
pumped into the z+ state, which is not being resonantly 
excited. When the initial state is the z+ electron state 
(dotted blue curve) the emission rate starts very low but 
rises slightly as a small amount of population is pumped 
into the z— electron state. Figure l^b) shows the aver¬ 
age number of photons detected D(T) as a function of 
the detection window duration for the two initial con¬ 
ditions. When the initial state is the z+ electron state 
(dotted blue curve), photon emission is relatively unlikely 
because the excitation laser is not resonant with the al¬ 
lowed (T_|_ transition; emission is still possible, however, 
due to weak remaining magnetic spin mixing. When the 


Time (ns) 



Detection window duration, T (ns) 

FIG. 5. (Color online) Spin read-out operation under a- 
resonant excitation. In (a) and (b) the two initial conditions 
are the z— electron state (solid red curve) and the z+ elec¬ 
tron state (dashed blue curve). (a) Photon emission rate after 
CT_ resonant excitation begins, (b) Average number of pho¬ 
tons detected as a function of detection window duration, (c) 
Fidelity of the spin measurement as a function of detection 
window duration. The vertical dashed line indicates the op¬ 
timum detection window. 


initial state is the z— electron state (solid red curve), 
photon emission is relatively likely. The fidelity of the 
spin measurement is A = (1 —pz+ +Pz-)/2 where Pz± is 
the probability of detecting at least one photon when the 
initial state is z±. The fidelity is plotted in Fig. He) as 
a function of detection window duration. The maximum 
fidelity occurs for a detection window of 165 ns, indi¬ 
cated by a vertical dashed line, beyond which the value 
of Pz- saturates but pz+ keeps increasing. For the op¬ 
timum detection window duration the calculated fidelity 
is 76.2%, which is slightly lower than the 82.3% mea¬ 
sured in the true Faraday configuration in Ref. [TH]. The 
fidelity would be improved with larger AG Stark shift, 
smaller Voigt magnetic field, smaller electron spin de¬ 
phasing and decay rates, or smaller trion dephasing rate. 
The values chosen here are either typical or feasible for 








11 


real experiments. 
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When the average number of photons detected is 
greater than unity, we can say that a single-shot mea¬ 
surement of the electron spin-state is possibl^i^. For 
these calculations, we have chosen an overall detection ef¬ 
ficiency e = 2.5% that results in an average number of de¬ 
tected photons greater than 1 for the optimum detection 
window duration; for T = 165 ns, D(T) = 1.01. Such 
an efficiency is relatively high, but should be achievable 
in recently developed “photonic trumpet” waveguide^, 
which have a demonstrated first-lens collection efficiency 
of 75%. Because application of a non-resonant laser is 
necessary for the AC Stark effect, a sample design such 
as the trumpet that achieves high collection efficiency 
through non-resonant effects is necessary, which elimi¬ 
nates most microcavities from consideration. 


IV. CONCLUSION 


We have described a scheme to accomplish a single¬ 
shot read-out of the spin-state of an electron trapped in a 
quantum dot while maintaining the capability to perform 
arbitrary coherent manipulation of the spin-state. A 
pseudo-Faraday configuration is produced by application 
of a Voigt geometry magnetic field and a far red-detuned, 
circularly-polarized laser that causes a spin-dependent 
AC Stark shift. The spin read-out is accomplished in 
the pseudo-Faraday configuration via spin-selective flu¬ 
orescence from spin-preserving cycling transitions. For 
typical quantum dot parameters and feasible detection 
efficiency, the spin measurement can be accomplished 
faster than the state is disturbed by the back-action, 
resulting in a single-shot read-out. The fidelity of the 
read-out is limited by the remaining spin-state mixing in 
the pseudo-Faraday configuration that is caused by the 
Voigt geometry magnetic field. Because the laser that 
produces the AC Stark shift can be switched on and off 
rapidly compared to the spin lifetime, this scheme offers 
the possibility to perform coherent spin manipulation in 
the Voigt configuration and then single-shot spin read¬ 
out in the pseudo-Faraday configuration. These capabil¬ 
ities comprise all three necessary single-qubit operations 
and will allow the investigation of more complex control 
and manipulation sequences. 


The authors wish to acknowledge helpful discussions 
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Appendix: Branching ratio 


The branching ratio for the pseudo-Faraday configura¬ 
tion can be derived by considering the four-level charged 
QD system interacting with the quantized multi-mode 
electromagnetic field. Extending the derivation in ref¬ 
erence [2H] to four levels, the state of the system as a 
function of time is 


|5»(t)) =C^{t) 1-04, vac) -h Csit) 1-03, vac) 

+ {CakW 102, k) + CikW 101, k)} (A.l) 

k 

and the interaction Hamiltonian in the rotating wave ap¬ 
proximation is 

V = + H.c.} (A.2) 

k i = l,2 

3=3,4 

where the |0i) are the eigenstates of the system, ciy = 
100 (0j| is the QD lowering operator, Ok is the photon 
annihilation operator for the electromagnetic field mode 
with wavevector k and frequency ojk, iVji is the transition 
frequency between states |0j) and \ipi), ftyk is a coupling 
constant, and H.c. means the Hermitian conjugate. The 
coupling constant is given by 


gij,k =-Vij ■ (A.3) 

where Vij = e {ipil r |0j) is the electric dipole matrix ele¬ 
ment in the eigenbasis and in general is a complex vector, 
Ck is the polarization vector of the k-mode, and 


fk = 



(A.4) 


where eg is the permittivity of free space, and V is the 
quantization volume. 

To determine the time evolution of the amplitude coef¬ 
ficients of |4>(t)) we substitute Eqn. (A.l) into the inter¬ 
action picture Schrodinger equation and equate the coef¬ 
ficients of similar kets. The result is a series of coupled 
linear differential equations: 









12 

Ca = 

(514,kC’ 

k 



.^.)t^2k) 

(A.5) 

C3 = 

(nts.ke’ 

k 



.^.)i^2k) 

(A.6) 

C'2k = 

— * ^2 (524,kC’ 

k 


+ ff23,ke-*(^^^-^'')‘C3) 

(A.7) 

Cik = 

— i ^2 (514,kC 


+ ffi3,ke-*(^" 


(A.8) 


These equations can be solved following the usual modes are closely spaced, which changes the sum over 
Weisskopf-Wigner theorjEH. First, we directly inte- k into an integral over k-space. The result is two cou- 


grate Eqns. ( A.7[ ) and ( |A.8| ) and substitute them into pled linear differential equations for C4 and C 3 , one of 

Eqns. (A.51 and ( |A.6[ ). Then we assume that the k- which is 




dt' 


Via ■ Ck 


'P 2 A ■ Ck 




(A.9) 


Each term of the integrand in Eqn. (A.9) has a factor 
of the form 


(u • 4) (v • ek) 


(A.IO) 


where e.g., u = V^a and v = 7^13. Since ek is a unit 
vector these factors only depend on the angular part of 
the integral over k-space. We define a general integral as 
follows: 


A = 


y d^f^k (u • 4) (v • Ck) (A.ll) 


which can be directly integrated to give 

, 47r 

A = — u • V 
3 


(A.12) 


Thus, the angular parts of the k-space integral in 
Eqn. (A.9) can be replaced using Eqn. (A.12) with ap¬ 


propriate substitutions. We can continue the Weisskopf- 
Wigner theory with the approximation that the coef¬ 
ficients CA(t) and C 3 (t) evolve much slower than the 
oscillation frequency oj (see reference [28] for details). 
This allows us to perform the integrals over k and f in 
Eqn. (A.9) to obtain 


= (r4i + r42) C4(t) - i (/331 + P 32 ) e*‘"«‘C3(t) 

(A.13) 


where the spontaneous decay rates r4i and r42 are 

(A.14) 


r ■■ = 


1 44 


V,. 


47reo3&3 

and the transition rates /331 and ,^32 are 

1 


Pji — 


47reo 37ic3 


4a;3- / ^ ^ 

'3 [^*1 ■ 


(A.15) 


where I = 4 wh en 7 = 3 and vice versa. An expression 
similar to Eqn. (A.13) can be obtained for C 3 by switch¬ 


ing all the indices 3 and 4. 

The transition rates Pji are zero when the dipole mo¬ 
ments for the two transitions to one electron state are 
orthogonal. This is the case for a charged self-assembled 
QD in the Earaday, pseudo-Faraday, or V oigt c onfigura- 
tions, which can be confirmed using Eqn. (A.20) and the 
eigenvectors of Hq from Eqn. ([^. Thus Eqn. (A.13) con¬ 
tains only the spontaneous decay rates. When the dipole 
moments are non-orthogonal, the re is t he possibility of 
spontaneously generated coherencd^^*^ and/or quantum 
interference caus ing p opulation transfer from one trion 
state to the othei®^^. 

The branching ratio is the ratio of the spontaneous 
decay rates F14 and r24. The difference between the 
transition frequencies ujai and a;42 is much smaller than 
their magnitude, so we can approximate the branching 
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ratio tb &s 


tb = 



(A.16) 


From the definition of Vij, above, we can use the com¬ 
pleteness of the Faraday basis {|a)} from Eqns. ([^ to 
obtain 


(A.17) 

a/3 


qi3j Qsij q24, and q42. We know that in the Faraday con¬ 
figuration the allowed transitions are circularly polarized 
when the photon is emitted in the z-direction. The dipole 
moments for such transitions are 


qi3 = ^ (^ + iy) 
q24 = iy) 


(A.18) 


Combining Eqns. (A.17), (A.18), and (i) we can define 
the electric dipole operator in the Faraday basis as 


d = qi3cr+ + q240-- + H.c. (A.19) 


where = e (q;| r |/3) is the electric dipole vector ma¬ 
trix element in the Faraday basis. Note that q^/j = qj^ 
and some of the q^,^ are zero. For example, q^a = 0 by 
parity for all values of a. Also, q43 = q2i = 0 at opti¬ 
cal frequencies. We know that the polarization selection 
rules of the Faraday conhguration are such that the al¬ 
lowed transitions are |1) O |3) and |2) o |4). Thus, we 
can assert that only 4 of the remaining q^^ are non-zero: 


and the transition dipole moments between the eigen¬ 
states are 


Vrj = (V-il d \i^j) 


(A.20) 


Thus, the branching ratio in Eqn. (12) can be derived by 
combining Eqns. (A.16) and (A.20). 
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